Angle-resolved photoemission spectroscopy of band tails in lightly doped cuprates 
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We amend ab initio strongly- correlated band structures by taking into account the band-tailing 
phenomenon in doped charge-transfer Mott-Hubbard insulators. We show that the photoemission 
from band tails accounts for sharp " quasi-particle" peaks, rapid loss of their intensities in some di- 
rections of the Brillouin zone ("Fermi- arcs") and high-energy "waterfall" anomalies as a consequence 
of matrix-element effects of disorder-localised states in the charge-transfer gap of doped cuprates. 
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Since the discovery of high-T c superconductivity 
in cuprates, angle-resolved photoemission spectroscopy 
(ARPES) has offered a tremendous advance into the 
understanding of their electronic structure How- 
ever, even though ARPES is continually strengthening 
our insights into the band structure and correlations 
in cuprates, it has also revealed many poorly under- 
stood phenomena, such as the incoherent "background", 
the sharp "quasi-particle" peaks near some points of 
the Brillouin zone (BZ), which form "arcs" of "Fermi 
surface" (FS) ([2] and references therein), widely stud- 
ied low-energy dip-hump and kink features (for review 
see [l|) and the more recently discovered steep down- 
turn of the dispersion toward higher energies (the so- 
called "waterfall") @, i, i, i, SHU- These anomalies 
have received quite different interpretations, involving, 
for exam ple, uncorrelated pj| and strongly-correlated 
@,[ll,[il,[ilEii.attice polarons, Migdal-Eliashberg-like 
approaches [3, Qj} , spinons and holons Q , spin polarons 
pj] , spin fluctuations [l(J EH and band-structure matrix 
element effects [l,[l9|. 

Many ARPES interpretations suggest a large FS (as 
an exception see e.g. |ll|) with nodal gapless quasiparti- 
cles, which are gapped or strongly damped in the antin- 
odal directions ((0,0) — > (tt, 0)) of the two-dimensional 
(2D) BZ. Importantly, extensive simulations of ARPES 
using the first-principles (LDA) band theory with the ma- 
trix elements properly taken into account [19[ reproduces 
well the topological features of momentum distribution 
curves (MDC), pointing to the large FS in optimally 
doped cuprates. However, LDA predicts that the un- 
doped parent cuprates are metallic with roughly the same 
large FS, while they are actually charge-transfer Mott- 
Hubbard insulators with the optical gap at 2 eV. This 
fact led to several powerful extensions of LDA, in par- 
ticular to LDA+U, which combines LDA eigenfunctions 
with strong Coulomb correlations introduced as a model 
parameter (Hubbard U) [20[, and the LDA+generalized 
tight-binding (GTB) method combining the exact diag- 
onalization of the intracell part of the Hamiltonian with 
relevant LDA eigenfunctions and Coulomb correlations 
and the perturbation treatment of the intercell hoppings 
and interactions [2~0 |. LDA-GTB Hamiltonian is reduced 



to the simpler effective t — J or t — J* model (t — J 
model plus three-center correlated hoppings |2l|) in the 
low-energy domain. 

LDA+GTB band structure of undoped cuprates with 
ab initio sets of tight-binding parameters [2l| describes 
remarkably well the optical gap, E ct « 2 eV both in 
antiferromagnetic and paramagnetic states of the un- 
doped La2Cu04. The valence band consists of a set 
of very narrow (< 1 eV) subbands where the highest 
one is dominated by the oxygen p states with the max- 
imum at k = g = (7r/2a, 7r/2a) (see Fig.l), while the 
bottom of the empty conduction band formed by d x i_ y i 
states of copper is found at (ir/a, 0). These locations of 
valence-band maximum and conduction-band minimum 
perfectly agree with ARPES intensity locus in hole doped 
La2-; C Sr ;E Cu04 and electron-doped Nd2-^Ce :E Cu04, re- 
spectively [22|. Importantly, the LDA+GTB approach 
predicts the charge-transfer gap at any doping with the 
chemical potential pinned near the top of the valence 
band (in hole doped cuprates) and near the bottom of 
the conduction band (in electron-doped cuprates) due to 
spin-polaron in — gap states. 

ARPES of undoped cuprates 0, i, @, S [H, [1| proved 
to be critical in the assessment of different theoretical 
approaches. It revealed an apparent contradiction with 
the t — J model. There is no sharp peak predicted by 
the model in undoped cuprates, but a slightly dispersive 
broad incoherent background, Fig. 2 (inset). Small lat- 
tice polarons due to a strong electron-phonon interaction 
(EPI) have been advocated as a plausible explanation of 
the discrepancy [13[. When EPI is strong, the spectral 
weight, Z, of the coherent small-polaron peak is very 
small, Z « 1 and, hence the peak can not be seen in 
experiment since all weight of the sharp resonance in the 
t — J model is transformed at strong EPI into the broad 
continuum. 

Unfortunately the energy distribution curves (EDC)in 
La2Cu04 , Fig. 2a, has only little if any resemblance 
to the small-polaron spectral function, which is roughly 
gaussian-like. Only by subtracting a "background" given 
by the spectrum near (ir/a, ir/a), Fig. 2a, one can account 
for the remaining EDC with the polaronic spectral func- 
tion [15]. This background problem obscures any reli- 
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FIG. 1: LDA+GTB valence band dispersion [2l[ amended 
with band tails (ladder lines) near T, (7r/2,7r/2) and (jr, tt) 
maxima ( here k is measured in 1/a) 



able interpretation of the broad ARPES intensities, espe- 
cially in underdoped cuprates, where the charge-transfer 
gap at 2 eV makes inelastic scattering events implausi- 
ble as an explanation of the background. Sharp peaks 
at (7r/2a, ir/2a) near the Fermi level, Fig. 2b, in doped 
cuprates also remains a puzzle. Small heavy polarons 
cannot screen EPI in lightly doped cuprates. Hence, if 
Z is small in the parent cuprate, it should also remain 
small at finite doping, so that the emergency of the peaks 
cannot be explained by a substantial increase of Z with 
doping. 

Here we show that amending the LDA+GTB band 
structure of doped cuprates by inevitable impurity band- 
tails, the ARPES puzzles as mentioned above are ex- 
plained. 

Doping of cuprates inserts a large number of impuri- 
ties into the parent lattice. Each impurity ion locally 
introduces a distinct level, E{, in the charge-transfer 
gap. The fact that the impurities are randomly dis- 
tributed in space causes the density of states (DOS) to 
tail, like in heavily doped semiconductors [25|. When 
there are many impurities within the range & of a lo- 
calised wave function ^(r), the random potential pro- 
duces low-energy states near maxima of the valence band 
at hole doping, Fig.l, or near minima of the conduction 
band at electron doping. As a result, ARPES intensity, 
/(k, E) = /b(k, E) -\-Ii m (k,E) comprises the band-tail 
intensity, Ij m (k, E) , due to localised states within the 
charge-transfer gap, and the valence band contribution, 
ib(k, E), of itinerant Bloch-like states. According to LDA 
band structures pj| the itinerant states are anisotropic- 
3D (specifically in La2Cu04) dispersing with c-axis k z 
over a few hundred meV. We suggest that this dispersion 
shapes the background making it so different from the 
incoherent background caused by EPI and/or spin fluc- 



tuations since k z is not conserved in ARPES experiments. 
On the other hand the incoherent background can be well 
described by a simple polaronic Gaussian in presumably 
more anisotropic insulating Ca2Cu02Cl2[24]. 

Here we focus on the band-tailing contribution de- 
scribed by the Fermi-Dirac golden rule as 

I im (k,E) = ^n(£)^|(^|A - V\iPi)f 5(E+A-Ei) 

Tin ^ 

(i) 

We define all energies relative to the chemical potential, 
/i, which is situated within the impurity band as shown in 
Fig.l. Only the impurity states with the binding energy 
E{ below \i = contribute at zero temperature. Here Ao 
is the amplitude of X-ray vector potential, and h—ks — 
1. 

We take the impurity wavefunction as [26], ^(r) = 
Fi(r)ip s (r), and the final state to be the normalised plane 
wave, ipf(r) — (Nv)~ x l 2 exp(ik • r). Here V ; g( r ) 1S the 
itinerant state at the top of the valence band, Fi(r) is a 
slowly varying envelope function, and N is the number 
of unit cells of volume v in the crystal. In the framework 
of GTB [21] one can expand ip s (r) using the Wannier 
orbitals, ip s (r) = N -1 / 2 J2 m w ( r ~ m ) exp(zg • m), and 
calculate the dipole matrix element in Eq.(l) as 

I im (k, E) = In(E) l/i(k - E)\ 2 S(E + A - E t ), (2) 

i 

where / = 27r(ed/m e ) 2 (Ao • k) 2 /v is propor- 
tional to the valence band matrix element squared, 
which is roughly a constant in a wide range of 
k near g, d = f drw (r) exp(ig • r), and /i(q) = 
(Nv) -1 J dr exp(iq • r)i^(r) is the Fourier transform of 
the impurity envelope function. 

Since the size of the envelope is large compared with 
the lattice constant, its Fourier transform strongly de- 
pends on q, which explains the experimental EDC 
and MDC as we show in the rest of the paper. 
We choose the impurity state to be hydrogen- like, 
Fi(r) = (Nv/ir^f) 1 / 2 exp(— r/&) as the hydrogen model 
accurately predicts many properties of shallow lev- 
els in heavily doped semiconductors, so that /i(q) = 
87r(^ 3 /7rA^v) 1 / 2 (l + 4 2 <^ 2 )" 2 for 3D impurity states, and 
fi(q) ex (1 + <? 2 £ 2 ) -3 / 2 for 2D states like localised surface 
states. It is important to recognise here that & is related 
to the impurity binding energy as £~ 2 = mEi, where m 
is roughly the hole effective mass. As a result we get 
J im (k, E) = xIn(E)M(k - g, E) with 



M(k-g,£) 



64tt 



(25 + A) 



5/2 



urn 3 / 2 [E + A+ (k-g) 2 /m] 



2/™i4^(£+A). 



(3) 



Here pim(E) = TV" 1 ^2 i 5(E — E i ) is the band-tail density 
of states (DOS) normalised to unity, and x = Ni/N is the 
impurity concentration per cell proportional to doping. 
In the 2D case the result is similar, M2£>(k, E) oc (E + 
A) 2 [E + A + (k,| - g f/m]- 3 p im (E + A). 
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FIG. 2: Band-tail EDC, Eq.(4), (solid lines) with pseudogap 
A = 300 meV and band-tail width 7 = 300 meV compared 
with relative EDC (symbols) near (it /2a, 7v/2a). Relative in- 
tensities are obtained by subtracting ARPES intensities of the 
parent compound, La2Cu04 (a), shifted by 5p, from EDC of 
slightly doped Lai.97Sro.o3Cu04 (b) as measured by Yoshida 
et al 2]. Both intensities have been normalised by their 
values at E — —800 meV and the chemical potential shift 
between two samples has been taken as 5/i — 70 meV. 



We notice that due to a very sharp dependence on q of 
the matrix element in Eq.(2) any uncertainty of k z does 
not smear out the strong dependence of Ii m (k, E) on the 
in-plane momentum component, kn. Averaging over k z 
simply replaces M(k — g, E) in Eq.(3) by 



M(k||-g,S) 



32c 



(E + A) 5 / 2 



vm [E + A + (kii - g; 



2/ m ]7/2 



p im (E+A), 
( 4 J 

where c is the c-axis lattice constant. Also M and M 
can be very large for shallow impurity states, M, M ^> 
1/x. Hence even the strong polaronic reduction of their 
weight, Z <C 1, does not make band-tails invisible in 
ARPES at finite doping, in contrast to a complete re- 
duction of the coherent band peak. 

Since the chemical potential shifts towards the band 
edge with doping, A in Eqs.(3,4) becomes smaller. 
Hence, the band-tail peak, Ji m (k, E), which is propor- 
tional to x, not only increases but also becomes sharper 
with doping as observed [2[ . To provide more insight into 
the shape and momentum dependence of experimental 
EDC we approximate the band-tail DOS by the simple 
form, p zrn (E) = [n/r(p/n+ l/n)](£/ 7 )^ exp(-£ n / 7 n ), 
where T(x) is the gamma-function. Exponents n,p de- 
pend on the dimensionality and the correlation length of 
the disorder potential: n = 2 both in 2D and 3D, p = 2 in 
2D and p = 7/2 in 3D for the long range random poten- 
tial correlations. In the short-range Gaussian- white-noise 
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FIG. 3: Waterfall effect in the band-tail ARPES intensity 
(white colour corresponds to the highest intensity). 
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FIG. 4: Real space Fourier transform (lower panel) of the 
square root of ARPES intensities (arb. units) at the Fermi 
level in Ca2-ccNa x Cu02Cl2 (upper panel, measured by Shen 
et al [23] for x = 0.12) reveals the real-space size (in units of 
a) of localised in-gap states. 



limit one obtains n = 1, 1/2 in 2D and 3D, respectively, 
and p = 3/2 in both dimensions [2jj| . We can separate im- 
purity and band contributions by subtracting normalised 
ARPES intensity of the parent cuprate from the intensity 
of the doped one. Then, the band-tail ARPES, Eq.(4), 
fits well with the experimental relative intensities at all 
momenta around g with m = m e , n = 2, and p = 7/2, 
Fig. 2. It describes the substantial loss of intensity with 
changing the momentum by only a few percent relative 
to g, as well as the shape of the relative EDC. 



4 



We argue that band-tailing can also contribute to the 
waterfall effect. There are impurity tails near local max- 
ima of the LDA+GTB valence band at T point (0,0) 
and at gi = (tr/a, 7r/a), as shown in Fig.l. Different 
from in-gap impurity states at g = (n/2a, it /2a) these 
localised states are hybridised with the valence band 
states of the same energy (shown by stars in Fig.l). 
However, the hybridisation could be insignificant, if the 
corresponding matrix elements of the random potential 
are small due to a large momentum separation between 
those states of the order of ir/2a. Hence, the impurity 
peaks reappear and disperse towards (0, 0) and gi at 
high bindingenergies, as observed in a number of doped 
cuprates 0,11, @, [E IS • We illustrate the waterfall in 
Fig. 3 by adding all three tail contributions, Ij m (k, E) oc 
n(E)[M(k^ £+£ 2 )+M(k|| -g, £7)+M(k|| - gl , E+E 2 )} 
where E 2 is roughly the valence band-width (we chose 
E 2 = 500 meV). We notice that the Fermi-Dirac distribu- 
tion, n(E), is replaced by its convolution with the Gaus- 
sian energy resolution function, n(E) — > [l — erf(E/S)]/2 
in plotting Figs. 2, 3 since the energy resolution 6 = 20 
meV is much larger than T « 2 meV. Also the photoe- 
mission intensity comprises both band-tail and valence 
band contributions, so that the resulting dispersion could 
be different from the anomalous band-tail dispersion of 
relative intensities, Fig. 2. 

Our theory proposes that the ARPES intensity near 
(7r/2a, 7r/2a) is proportional to the square of the Fourier 
component, /i(q), of the impurity wave- function enve- 
lope, Eq.(2). Therefore, we can find the real-space im- 
age of the function, Fi(r), by taking the Fourier trans- 
form of the square root of the experimental intensi- 
ties, Fig. 4 (upper panel). Here we show the intensities 



near the Fermi level measured in Ca2-^Na :E Cu O2CI2 
[231 ]. which are very similar, if not identical to those in 
La 2 -xSr a; Cu04 (compare Fig.l [23] and Fig.2 in [2]). The 
real-space image (lower panel, Fig. 4) reveals some band- 
mass anisotropy and the size of the localised state of 
about 20 lattice constants justifies the "envelope" ap- 
proximation [26[ used for the impurity wavefunction. 

In summary, we have proposed an explanation for 
sharp " quasi-particle" peaks, "Fermi-arcs", and the high- 
energy waterfall in cuprates as a consequence of matrix- 
element effects of disorder-localised band-tails in the 
charge-transfer gap of doped Mott-Hubbard insulators. 
Importantly if holes are bound into bipolarons, the chem- 
ical potential remains within the single-particle band-tail 
at the bipolaron mobility edge even up to optimum dop- 
ing, in agreement with S — N — S tunnelling experiments 
[28[ and insulating-like low-temperature resistivity of un- 
derdoped cuprates. In this case A in Fig.l is half of the 
bipolaron binding energy [11], which is also the normal 
state pseudogap [29]. Recent scanning tunnelling mi- 
croscopy at the atomic scale found intense nanoscale dis- 
order in high-Tc superconductor Bi 2 Sr 2 CaCu208+<5 [!o[ 
telling us that band-tailing indeed plays the important 
role in shaping single-particle spectral functions of doped 
Mott insulators. 
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